Cancellation of atmospheric turbulence angle-of-arrival fluctuations 



Marcelo V. da Cunha PereiraQ Lufsa A. P. Filpi, and C. H. Monken 

Departamento de Fisica, Universidade Federal de Minas Germs, 
Caixa Postal 702, Belo Homzonte, MG 30123-970, Brazil 
(Dated: February 16, 2012) 

By using a turbulence chamber which emulates the dynamic behavior of the actual turbulent 
atmosphere, we show that the most deleterious wavefront aberration imparted to an optical beam, 
the angle-of-arrival fluctuation, which corresponds in practical situations to more than 87% of all 
the wavefront aberrations in free-space propagation, can be cancelled using two-photon correlation 
measurements and passive linear optics. 

PACS numbers: 42.68 Ay, 42.50 Ex, 03.67 Hk, 42.25 Dd 



I. INTRODUCTION 

Usually, when light is sent through free-space links, 
information is encoded in polarization because clear-air 
propagation disturbances don't affect it much [1] but, 
because it is a two-dimensional degree of freedom, the 
amount of information which can be transmitted per pho- 
ton is limited. Photons spatial degrees of freedom allows 
for much larger dimensionality space to be accessed [SHE] , 
the problem being that information encoded in these are 
typically very degraded by the atmosphere [7] . 

As has been shown recently in [8], the strong correla- 
tions in entangled states can assist in performing clas- 
sical communication under the effect of specific sources 
of error. In this letter we report an open possibility of 
encoding information in the spatial degrees of freedom of 
a two-photon beam in a way to make it resilient against 
atmospheric turbulence or, more specifically, wavefront 
tilt caused by atmospheric inhomogeneities. Among the 
different phase aberrations imparted by the refractive in- 
dex fluctuations in turbulent atmosphere, this has long 
been known [9] for being the most deleterious in terms 
of loss in resolution when long-exposure imaging is used 
[TU] , accounting for as much as 87% of the total wavefront 
distortion. 

This aberration is routinely compensated by means of 
adaptive optics, where actuators change the shape of a 
telescope depending on a probe feedback [TT] . The differ- 
ence between this and our approach is that ours use the 
inherent correlation of the twin photons in such a way 
that no active compensation is required. We also don't 
require a bright source as in schemes that employ optical 
phase conjugation |12j . 

II. THEORY 

We consider a situation in which a beam of wavelength 
propagates in the z direction for a distance D, over which 
it suffers the effects of random refractive index fluctu- 
ations, and then is focused by a positive lens of focal 



length / onto a photodetector. The beam may come 
from a laser, in which case we'll be interested in its in- 
tensity at the lens focal plane, or from a spontaneous 
parametric down-conversion (SPDC) process, where a 
nonlinear crystal is pumped by a well collimated laser 
beam of wavevector k p . For the latter situation, we'll 
measure the fourth-order correlation by using coincidence 
measurements and the entity that behaves as a beam 
is the correlation of the photon pairs itself. The key 
result in this letter is to show that when dealing with 
SPDC, a proper preparation of the down-converted pho- 
tons, through a passive and linear setup, can mitigate 
the effects of wavefront tilt. Transverse coordinates in 
the lens plane are denoted by p while those in the plane 
of detection are u, with eventual subindexes 1, 2 labeling 
different eletromagnetic modes. 

The refractive index fluctuations can induce both am- 
plitude and phase fluctuations in the beam, the latter be- 
ing the one in which we are interested here; in particular, 
wavefront tilt. It is described by an extra multiplicative 
factor of exp \ik6j ■ Pj\ in each of the field mode opera- 
tors, where 6 is a random variable denoting the tilt angle. 
This leads to the following two-photon wavefunction in 
the plane of the lens [5J: 

A( Px , P2 )=e(^±^) V( Pl - P2 ) 

X e <TF«iPl e <¥fl»fti , (i) 

where E(p) is the transverse profile of the pump in the 
lens plane in the absence of turbulence and V(p) is 
the propagated phase-matching function, whose Fourier 
transform will be given below. Now, in the focal plane 
of the lens, the field amplitude is the Fourier transform 
of A(p 1 ,p 2 ), from variables pj to q., = |jUj (j = 1,2). 
The linearity of the extra phase factor in p means that 
the resulting transform can be put in terms of only the 
variables q+ = qi + q2 and q = qi ~ q2 : 

5+ ( qi , q 2 ) cx E (q+ + k p [6 1 + 9 2 ]) V (q- + k p - ~ - j 
^E(q + + k. p [e 1 +e 2 ]) . (2) 
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In the above expression, V(Q) is the function 
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sine ( , L I Q 1 2 ) propagated from the crystal's output 

\ Kp Tip J 

face up to the lens plane and also the Fourier transform 
of V up to some phase factor [5]. In the thin crystal 
regime, it varies much slower with its argument than E 
and for this reason it was considered to be a constant. 

If we get one of the photons of a pair (say, photon 2) to 
suffer an inversion in transverse coordinates p 2 — > —P2i 
the two-photon amplitude function in the focal plane of 
the lens will be, instead, 



g-(qi, 02) oc E ( q + + k p 1 2 2 j . (3) 

Whatever the case, the coincidences rate is given by 
R c = ( \ g \ 2 ), with the average being taken over the tur- 
bulent air variables. It can be seen that if 0\ and 62 are 
correlated, cancellation will occur in g_ (while the tur- 
bulence effects on g + will be the same as those on a laser 
beam with wavevector 2k p ). To enforce this, entangled 
photons need to be transversally confined to an area of 
the order of g 2 , where g c is turbulence coherence length 
and gives the distance above which phase aberrations be- 
come uncorrelated. Of course this has to be mantained 
during the entire propagation, since the net distortion is, 
as opposed to what might be suggested by the adopted 
model, a result of accumulated small perturbations over 
the whole traversed length. 



III. EXPERIMENT 

To emulate the conditions present in the actual atmo- 
sphere we adapt a model put forward by Keskin |13j and 
depicted in figure [T] it is an aluminum box inside which 
air of different temperatures are mixed, which results in 
a random temperature field leading to an index of refrac- 
tion that varies randomly in space and time |14j . Air is 
blown into the chamber by two fans, one of which has 
a set of resistors (6.4 f2) in front of it. The air going 
through the resistors gets heated up and proceeds to the 
center of the chamber where it mixes with the air blown 
by the other fan, which is at room temperature. The op- 
tical beam goes in and out of the chamber through two 
oppositely positioned holes, propagating perpendicularly 
to the direction of the fans and through the center where 
mixing takes place. Two other fans work as exhausts to 
the mixed air in order to ensure the stationarity of the 
process. 

Our source is a 5 mm long BiBO crystal cut for type-I 
(it generates photons with same polarizations) collincar 
SPDC (NL Crystal in figure |2|. It is pumped by a cw 
He-Cd laser (wavelength 325 nm, power 42 mW), gen- 
erating entangled photons the polarizations of which are 
rotated by a half-wave plate to 45° relative to the hori- 
zontal plane. Some of the pairs will have their component 
photons leaving at different directions off the subsequent 
polarizing beam-splitter (PBS) and will go through dif- 
ferent arms of an interferometer before going into the 
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FIG. 1. (color online). A turbulence chamber that emulates 
the turbulent mixing of air at different temperatures which 
occurs in the atmosphere. 



turbulence chamber. Because of another PBS right be- 
fore the detectors, only these pairs will give rise to coin- 
cidence counts; pairs whose photons go through the same 
arm of the interferometer don't give rise to coincidences 
and are thus discarded. Inside the interferometer a half 
and a quarter-wave plates are placed oriented at 45° rel- 
ative to the horizontal plane, so that a photon that was 
initially reflected will be transmitted by the PBSs and 
vice-versa. This results in all sampled photon pairs leav- 
ing through the same port of the second PBS, but with 
one photon in each pair suffering an additional reflection, 
thus effecting the inversion pi — > —pi- After collinearly 
propagating through the turbulence chamber, they are 
collected by a lens of focal length 50 cm that focuses 
them on the detection planes. The detectors are Per kin 
Elmer SPCM-AQR-14 photon counting modules in front 
of which we place 10 nm bandpass filters centered at 
650 nm and 50 \wn wide slits. The slits, each having one 
of its dimensions much larger than the transverse widths 
of the beam, effectively average out the effects pertaining 
to the y direction, leaving as parsed only the x direction, 
along which the inversion scheme is implemented. 

In our experiment, a laser beam with a gaussian pro- 
file pumps the nonlinear crystal leading to an also gaus- 
sian coincidence profile (see equations ^ and (J3j) ) . We 
measure the curve with 20s for each point and perform 
a gaussian fit, repeating the process for different turbu- 
lence strengths. This can be varied by altering the in- 
put voltage in the resistors. The widths of the fitted 
curves are then registered and plotted against the square 
root of the Rytov variance \J a%, a conventional figure 
of merit for the turbulence strength. The measurement 
is repeated by removing the interferometer and placing 
a half-wave plate just before the last PBS, to ensure co- 
incidence counts. These two sets of measurements corre- 
spond to the g- and 5+ cases, respectively, and are shown 
together in figure [3] Another measurement is performed 
using a cw He-Ne laser beam propagating under the same 
conditions, and the collected data is also included in the 
same plot. One can see smaller broadening in the g_ case 
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following formula (T5] : 
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1.33A 5/6 



(4) 



, where A = ^5 and cto = cr (0)- At the highest power we 
get a Rytov variance of a 2 R = 0.17, close to the value 0.3 
that conventionally separates weak and strong turbulence 
regimes pQ. It is an established fact that the theory de- 
scribing phase fluctuations in turbulence has its validity 
extending far into the strong regime, so that the cancel- 
lation exposed here should be expected also under this 
circumstance fl4l. 



FIG. 2. (color online). Schematics of the experimental setup. 
The interferometer following the entangled photons source 
(NL Crystal) ensures that one of the photons suffers an addi- 
tional reflection relative to the other. 
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FIG. 3. (color online). Gaussian widths of the measured pro- 
files as a function of 



' a\ (see text). Blue circles and red 
triangles are for coincidence measurements with and without 
the interferometer, and green empty squares is for the inten- 
sity measurements of a He-Ne laser beam. 



as well as similar broadening in g + and laser beam data, 
as predicted. 

It is, however, to be kept in mind that wavefront tilt 
is only but one of the phase aberrations imparted to the 
beam and that the next relevant ones, the quadratic aber- 
rations, are not cancelled in the present scheme, and in- 
deed start showing up at higher voltages. Their impor- 
tance are, however, overestimated in this experiment be- 
cause of the limited dimensions of the turbulence cham- 
ber. A more complete description of the turbulent air 
effects has to take into account its size and, as shown in 
[15j . a consequence of this confinement is a shift of the 
relative importance of the degrading effects from wave- 
front tilt towards those of higher-order aberrations. 

Calibration of the turbulence chamber was made by 
sending an expanded He-Ne laser beam through it, and 
focusing onto a CCD camera. Using the known beam 
radius w — 68/im and traversed path length inside the 
chamber D = 76 cm, and measuring the widths <r(P) at 
the focal plane as a function of power P dissipated by 
the resistors, the Rytov variance was computed with the 



IV. DISCUSSION 

The cancellation can be understood in terms of 
Klyshko "advanced wave" picture [T7], a formal equiv- 
alence between coincidence detection of down- converted 
photons and the detection of a photon coming from one 
of the detectors, reflecting at the non-linear crystal and 
propagating back to the other detector. In the latter sce- 
nario the cancellation can be seen as follows: the photon 
leaves the first detector, has its direction changed by the 
turbulence, inverts the deflection angle at the interfer- 
ometer and then get it undone upon propagating again 
through the same inhomogeneities inside the chamber be- 
fore reaching the second detector. 
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FIG. 4. (color online). Unfolded "Klyshko" picture of the can- 
cellation process. Photon has turbulence deflection inverted 
in the preparation setup and then undone when going again 
through the chamber. Solid line shows path with deflection 
by turbulence inhomogeneities and dashed line shows path 
without deflection. 

An interesting question arises as whereas entanglement 
is actually necessary to attain cancellation: in fact what 
was needed here was a u\ + 112 dependence on the de- 
tectors' positions u\ and U2, the kind of which appears 
on the conditional fringes of Greenberger et al. [IB]. In 
[19] it is shown that for a known transfer function, clas- 
sical correlations can be engineered in such a way as to 
reproduce any joint detection probability attainable in a 
single plane with entangled photons. Here, we only make 
use of the far-field correlations, but the transfer function 
is random and unknown, and it is not clear if in this case 
one can reproduce the obtained results with only classical 
correlations and a passive setup. 
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Sending unspoiled entangled photons, however, is an 
advantage in itself because it allows for the use of quan- 
tum communication protocols such as super-dense cod- 
ing and entanglement purification. Investigations on the 
behavior of spatial entanglement under atmospheric tur- 
bulence has been already pursued but the question still 
remains as to whereas one can effectively protect it by 
properly preparing the states to be sent. This is to be 
addressed in future research. 

It should also be stressed out that the tilt cancellation 
can also be verified theoretically by using the generalized 
Huygens-Fresnel method, which includes in the propaga- 
tors the effect of turbulence 1 20 . This method have been 



used somewhat extensively together with an approxima- 
tion usually only employed to make calculations more 
tractable, which is assuming a square scaling-law for the 
turbulence structure functions instead of the more correct 
Kolmogorov 5/3 law. As shown in |21j . it happens that 
this approximation means physically that only angle-of- 
arrival aberrations are taken into account. 
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